Time-dependent theory of non-linear response and current fluctuations 
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A general non-linear response theory is derived for an arbitrary time-dependent Hamiltonian, not 
necessarily obeying time-reversal symmetry. This allows us to obtain a greatly generalized Kubo 
type formula. Applied to a mesoscopic system with any type of interactions, and coupled to multiple 
probes and gates with arbitrarily time-dependent voltages, we derive current-conserving differential 
conductance and current fluctuation matrices obeying a generalized Fluctuation-Dissipation Theo- 
rem. This relation provides a common explanation for asymmetries of the excess noise in several 
non-linear mesoscopic systems, as well as of its surprising negative sign. 



The linear response theory is a cornerstone of the quan- 
tum theory. It allows deriving extremely useful formulas 
such as the Kubo formula or the Fluctuation- Dissipa- 
tion Theorem (FDT) which are extensively used in all 
fields. It has been somewhat a common belief that the 
validity of the Kubo formula is limited essentially to the 
linear regime, motivating other alternatives in correlated 
systems [ij. In the present work we show that, in fact, it 
is possible to greatly generalize the Kubo formula within 
a non-linear response theory capable of addressing non- 
equilibrium situations and to include interactions as well 
as any (possibly time-dependent) Hamiltonian. Our work 
goes beyond and is more general than previous extensions 
performed in the sixties 'T\ and more recently in the sta- 
tionary regime[3,]: these works were not adapted to a 
mesoscopic context and to time-dependent Hamiltonian 
or/and voltages, and even under these limitations we are 
able to yield a much more compact proof. 

Even though our theory is not restricted to condensed 
matter theory, we choose to apply it to develop a novel 
transport formalism for arbitrary mesoscopic systems 
connected to two or many probes with any time depen- 
dence of their electrochemical potentials. Such multi- 
probe geometries are of great interest for e.g. reveal- 
ing statistics and entanglement, as in Hanbury-Brown 
Twiss setups or Mach-Zhender interferometers. For these 
systems, we express the differential out-of-equilibrium 
conductance matrix in a microscopic way, ensuring the 
current conservation and gauge invariance. We thereby 
provide a convenient general framework to describe fi- 
nite frequency and/or time-dependent behavior of these 
multi-probe systems which was lacking up to now. We 
also solve subtleties which have been a subject of de- 
bates since the development of the scattering approach 
[3, @ • Our formalism is an important achievement not 
only for time-dependent voltages or/and Hamiltonian, 
but is already crucial for the stationary regime. It offers a 
promising alternative to other approaches [11, 0, H, 0, 0, [1] , 
which even though successful and extended to nonlinear- 
ities and/or AC voltages are not suited to deal in a 
systematic way with a majority of the strongly correlated 
systems [l0| . 

Since we provide for any time-dependence of voltages, 
we can consider not only AC voltages [ll| to study photo- 
assisted transport [H, [iIjH, but also injection of elec- 



trons on demand [15|, |l6|, |l7[ , classical sources of noise, 
pumping, or mixing setups where the potentials in the 
reservoirs or gates have different periods, etc. We can 
also consider spontaneous generation under a DC bias, 
such as finite frequency (FF) noise [j], or combine both, 
for instance by applying time-dependent voltages and 
consider FF current fluctuations, which in this situa- 
tion depend on two frequencies [18| and form a matrix 
containing both auto-correlations and cross-correlations. 
Using our out-of-equilibrium and time-dependent Kubo 
formula, we show that such a matrix obeys a general 
time-dependent out-of-equilibrium FDT. Then we dis- 
cuss its important consequences in the limit of station- 
ary Hamiltonian and voltages. A first application of the 
present formalism has been done in quantum wires and 
carbon nanotubes p^ . 

For generality, we consider a system with an ar- 
bitrary time-dependent Hamiltonian Ti(t), which in- 
cludes any interactions or disorder, and does not nec- 
essarily obey time- reversal symmetry. 'H{t) depends 
in linear/non- linear and in local/non-local way on a 
set of time-dependent parameters generically denoted 
by X{t'). We express the functional derivative of the 
average of any operator O at time t with respect to 
X{t'). For this purpose, the Hamiltonian is split into 
the part which does not depend on X, denoted by 
7io(i), and another which depends on X: Ti.{t) — 
7io(i) + 'Hi{t,X). One switches to the interaction pic- 
ture where Tii is viewed as the interaction Hamilto- 
nian. Then 0™*(t) = Uo{-oo, t)0{t)Uo{t, -oo) where 
iUdtUo{t, — oo) = 'Ho{t)Uo{t, — oo). Even though not nec- 
essary, we prefer to exploit the Keldysh formulation to 
make our argument more compact. The Keldysh time 
contour has two branches labeled by 77 = ±, going from 
—00 to 00 on the upper one and inversely on the lower 
one. Tk is the Keldysh ordering operators which makes 
time (anti-time) ordering on the upper (lower) contour, 
while operators labeled by — are on the left to these la- 
beled by +: {TKA+{t)B-{t')} = {B{t')A{t)). O can be 
labeled by r/ = -I- or — for its average expression, of which 
the functional derivative with respect to X{t') is: 
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we obtain the central result of this paper: 



6X{t') 



/ dre{t 



n 



0{t). 



5X{t') 



•(1) 



Thus the functional derivative of the average of an oper- 
ator O can be expressed in terms of its commutator with 
the generalized force, defined by the functional differen- 
tial of 7i(i). The average is taken in the presence of TL{t) 
and with an initial density matrix which has not to be 
thermal. We will express higher order differentials, and 
allow for O to depend on X separately [20j . 

Now we apply this formula to a mesoscopic system 
connected to N terminals with electro-chemical potential 
A*n(i) = sVn{t) and a total charge operator Qn, the cur- 
rent operator /„ for each n = 1, .., N being /„ = ~dtQn- 
The system is described by TCq {t) , which can for instance 
include time-dependent scatterers, to which we add the 
coupling to terminals [1, : 



M 



(2) 



It is possible to justify this coupling [20] by generaliz- 
ing the formalism developed in one-dimensional wires 
[Tg . ITgj . We include the gates in the N terminals. We 
do not use the currently adopted coupling in terms of 
capacitances. Actually, this would be too specific as 
the latter could acquire a frequency dependence as well 
[20j . It is now possible to express in a microscopic 
way the differential out-of-equilibrium conductance with- 
out any constraint on the Hamiltonian, the initial den- 
sity matrix. We use Eq.(IT]) where X{t) is replaced by 
Vn{t) on which the interaction Hamiltonian in Eq.Q 
depends now in a local and linear way. Denoting by 




FIG. 1; A mesoscopic system coupled to many terminals in- 
cluding gates with arbitrary time-dependence of their electro- 
chemical potentials. The time-dependent Hamiltonian Ti. de- 
scribes any interactions or disorder and can depend on other 
parameters X{t') in a non-local and nonlinear way. Differ- 
ential of the current average at terminal n either with 
respect to V„i (t') or to X{t') can be expressed through a gen- 
eralized response formula keeping all Vn and X finite. 



Gnu'{t,t') = S{In{t))/SVn'{t') for n,n' - l..iV, we get 
the crucial result: 



hGnn'{t,t') = e{t - t') {[I„it), Qn'it')]) , 



(3) 



where the average is calculated still in the presence of 
the time-dependent Hamiltonian Ho{t)+7ii{t). Let's de- 
note by F{lu,lu') ^ J J dtdt'e"^*-"^'*' F{t,t'), the double- 
Fourier transform of a function F. Then Grm' (w, w') is 
the variation of (/„(tt')) when an infinitesimal modulation 
at uj' is added to K'(i') ^ F„/(t')+'y„'(t^')e'"'*', keeping 
Vn' (f) finite. 

In particular, Eq.Q should apply to interacting quan- 
tum dots weakly coupled to many reservoirs. In this case, 
it was claimed that the Kubo formula was not appropri- 
ate, contrary to what we have shown here, and the so- 
called non-equilibrium Green's function formalism was 
proposed as an alternative [Ij . In order to ensure the cur- 
rent conservation and gauge invariance which it missed, 
it has been adapted to a multi probe geometry where the 
values of the currents were modified Q. We think that 
this procedure lacks transparency and generality. 

Here, by carrying a total charge opposite to that 
on the system, the gates ensure the conservation law: 
Qn = 0- This guarantees gauge invariance auto- 
matically: a translation of all potentials Vn{t) by the 
same function V{t) has no effect on Ti,i{t) in Eq.®. 
Using Eq.(I21), we get simultaneously the two constraints 
I]„ G„„'(tJ, w') = X]„' G'„„'(w, w') = 0, the second one 
corresponds to gauge invariance. Thus the latter is not 
required independently of the current conservation, con- 
trary to previous works in the scattering approach. We 
emphasize that there is not necessarily time-reversal sym- 
metry, thus no symmetry of the matrix G. It is not 
hermitian neither. One can show that : G^^,{uj,uj') — 
G„„'(— a;, — tj'), where the star stands for the complex 
conjugate. 

Now we consider the non-symmetrized current fluctu- 
ations matrix whose elements are given by: 

Snn'{t,t') = {In'{t')In{t)) " {In'{t')) (/„(*)) . (4) 



Let's consider: 



S±(i,t') S{t,t')±S^ {t',t), 



(5) 



the symmetric and the antisymmetric parts of the current 
fluctuation matrix, where the subscript T refers to the 
transpose. Notice that one has S'„'„(i',t) = S'^^,{t,t'). 
One can show easily that fidt'Gnn'{t,t') — hdtGn'n{t',t) = 
S~^,{t,t'), which, once Fourier transformed, gives (see 
Eqs.jlEl)): 

S"(a;,w') = ^^^'(a;,^') - S'„/„(-cj', -cj) 

= - WG(w, uj') - hLuG''{Lu', LU). (6) 

This is a novel FDT which extends to out-of-equilibrium, 
time-dependent Hamiltonian and voltages, and with or 
without time-reversal symmetry. It allows to relate 
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both symmetrized and non-symmetrized fluctuations (see 
Eq.([5])), as one can inject this expression into the r. h. 
s. of: 



2S = 



(7) 



One can show that Sn'n{uj' ,uj) — S'*„/(aj, w'). Let's spec- 
ify to w — Lo' , which amounts to integrate out over t—t' to 
get the DC component of the fluctuations in the presence 
of any set of Vn{t). Eq.® reduces to : 



S-{uj,uj) = -2hujG''{uj,uj), 



(8) 



where the Hermitian part of the matrix G is given by: 



2G'' = G + Gt. 



(9) 



We note that S(w, w) and hence S^{uj,uj) become hermi- 
tian, which is not the case for G(w, lo). Only the diagonal 
elements of these four matrices are real. 

Let us now comment briefly on the case of a periodic 
potential. Allowing for generality a different frequency 
n„ in each terminal n as in mixing setups, one requires 
that at least for one terminal n one has lo — lo' = Ifln 
with / an integer. But both uj and uj' can be different 
from a multiple of all ri„ (2l| . 

We will focus in the following on both time- 
independent Hamiltonian and potentials in the reservoirs 
[2^ . Then invariance by time translation is restored and 
one requires lo' = u in Eqs. pi6p . Let's keep similar no- 
tations for F = S, G but stress their dependence on the 
voltage vector V = {Vi,V2.-Vn): 

F(w',w) = (5(w' - w)Fv(w). 

We also introduce the "excess AC differential conduc- 
tance" and excess FF noise matrices for later use by: 



AFv{uj) = Fyr{uj) - Fv=o(t^) 



(10) 



As noticed above for the case lo = lu' , the matrices Sy, Sy 
become hermitian while the matrix Gy is not. The asym- 
metry of the FF fluctuation matrix is analogous to Eq. ([5]) 
with a unique frequency (see Eq.([9])): 



Sv(cj) = -2ncjG^(a 



(11) 



This generalizes the scalar noise asymmetry obtained 
in Refs.'^, [2^. It yields directly the standard equilib- 
rium FDT provided time-reversal symmetry is ensured. 
In this case the equilibrium noise matrix obeys the de- 
tailed balance equation: Sv=o{—(^) = e^'^Sv=o{i^)- Thus 
Eqs.iniSD yield: Sy=o(t^) = 2;itjiV(w)G^^o(w), where 
N{uj) ~ l/(— 1-f e'^'^). Equation offers a generalized 
FDT which extends to both out-of-equilibrium and/or 
violation of time-reversal symmetry. 

It has as well other useful consequences. In case 
Gy(w) is known, Eq.© relates the symmetrized to non- 
symmetrized current fluctuations. Similarly, the fluctua- 
tions for negative (respectively positive) frequencies can 
be deduced from those at positive (respectively negative) 



frequencies. A more interesting alternative is to deduce 
Gy(w) from S~(a;), which allows to get rid of any back- 
ground undesirable noise, being a difference, and is not 
subject to the limitations on frequencies as the AC con- 
ductance. The latter are due to capacitive effects and to 
the equilibration condition in the reservoirs: <C 1 

where Ti„ is the inelastic time^in order to deflne a quasi- 
equilibrium distribution [l^, [l^, . 

Another important feature which Eq. pT)) clarifies con- 
cerns the asymmetry of the excess FF noise, Ea.([TUl). 
While many theoretical works used to study the sym- 
metrized noise, it turns out that most experiments 
are based on quantum detectors measuring the non- 
symmetrized excess noise [1^ H^, which has been the 
subject of few theoretical works with correlation effects 
[3, [13, HI]. An important question is under which 
criteria one can violate the symmetry obtained in a ba- 
sic coherent conductor p|, thus giving an evidence for a 
quantum measurement. It is interesting to discuss the 
asymmetry of the full excess fluctuation matrix within 
our formalism, thus to find the criteria for that of excess 
cross-correlations as well. This can be achieved by using 
simply Eq.ini) (see Eqs. ([51111)) ) (29j: 

ASv(-w) - AS'^(w) = 2hLoAG'{,{uj). (12) 

The asymmetry between AS'„„'(— w) and A5'„'„(u;) re- 
quires that AGJJ^, ^ 0. This yields a necessary criteria: 
non-linearity! However this is not sufficient for different 
terminals n ^ n': one could have G^^, = at any V, 
thus get symmetry of excess cross-correlations even with 
non-linearity, which clarifies this fact obtained in chiral 
edges of the FQHE [2^. At the same time, excess auto- 
correlations were found to be asymmetric in this system, 
as well as in quantum wires and carbon nanotubes p^ . 
Interactions are however not necessary for that. For in- 
stance, asymmetry holds in noninteracting systems where 
the transmission is energy dependent, such as a wire 
with two barriers (soj , hybrid structures 
son junctions measured experimentally 
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and Joseph- 
Again non- 
linearity is a common origin for the asymmetry. 

A related fact can be shown in the case of a tunneling 
junction with arbitrary interactions or disorder. For lu — 
qV/fi ^ ksT^ q being the effective charge, we can show 
that ASnn{uj) = [20| while A5'„„(— w) = ?ia;AG„„(ijj) 
do not vanish in non-linear systems. This generalizes and 
explains such behavior obtained in Refs.[ig.[2^. 

Considering again any non-linear system, it is interest- 
ing to introduce the combination we call the "modified 
excess noise" : 

ASy(w) = Sv{lo) - 2hLoN{Lu)G'^{Lo). (13) 

It restores symmetry with respect to any posi- 
tive/negative frequencies: ASy(— w) = AS'y(ijj). In 
view of Eq.lHI]), AS'y(cj) = (1 + iV(w))Sy(w) - 
N{Lu)Svi—(^)- 1^ ^ linear system, this becomes the ex- 
cess noise as one can check using the standard equilibrium 
FDT, and is identical to the combination measured by a 
detector [23l] having the same temperature as the system. 
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Now we show how the out-equihbrium FDT, Eq. ([TT|) . 
solves the paradox of the negative sign of the excess noise, 
which looks counter-intuitive as applying a bias is ex- 
pected to induce more noise, thus the nomination "ex- 
cess" . We focus on auto-correlations as they become real 
and can be interpreted in terms of emission/absorption 
spectrum. In two-terminal geometries, they could be neg- 
ative, such as in Luttinger liquids whether symmetriza- 
tion is performed [H, [33| or not [l9|, [ll] , or without in- 
teractions for an energy-dependent transmission [sil, [13] ■ 
Indeed, for huj ^ kgT, the equilibrium noise vanishes 
thus AS'„„(ijj) — S'„„(w) 0], which can be shown, by a 
spectral decomposition, to be always positive being the 
correlator of the same current at terminal n. But the 
absorption excess noise can be negative if AG'Jj„(a;) is 
negative enough, see Eq. (jl2p . Symmetrized excess noise 
contains both emission and absorption, thus can be neg- 
ative too. 

To conclude, we have derived a general time-dependent 
response formula for an arbitrary Hamiltonian depending 
in a local/non-local and linear/nonlinear way on time- 
dependent parameters. This yields a microscopic and 
current-conserving expression of the differential conduc- 
tance matrix in a multi-probe mesoscopic system with 



arbitrary time-dependence of the Hamiltonian including 
the electrochemical potentials in the terminals. We have 
deduced a general time-dependent out-of-equilibrium 
FDT which yields in particular the extension of the equi- 
librium FDT in case time-reversal symmetry is broken. 
Its application to the stationary regime has shed light on 
the asymmetry and sign of the excess FF noise in non- 
linear systems. 

Besides these applications, and that operated in quan- 
tum wires [l9l | , our theory offers a new promising frame- 
work to study systematically time-dependent transport 
in non-linear systems, including strongly correlated ones, 
and to consider pumping or mixing setups for instance. 
Remarkably, we are able to obtain as well higher or- 
der differential of any {0{t)) with respect to many time- 
dependent parameters, and therefore that of the current 
in probe ni with respect to the potentials in n2..nM, as 
will be reported elsewhere ^2^. 
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